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TECHNICAL MEMORANDUM X-53984 


CALCULATIONS OF ASTROPHYS ICAL 
PARTITION FUNCTIONS 

SUMMARY 

This report outlines a method for the calculation of the internal energy 
partition function of an atom or ion immersed in a plasma. The influence of 
the surrounding plasma on the atom or ion is considered by including a 
probability function P. in the sum over the energy states: 

00 -X./kT 

partition function U = ^ g. P. e 1 

i=i 1 1 

In this expression g. is the statistical weight of the ith level with excitation 

energy X-> T is the absolute temperature, and k is the Boltzmann constant. 

The function P^ denotes the probability that an atom or ion in the ith energy 

level is not perturbed by surrounding ions, and it is a function of the tempera- 
ture and electron pressure of the plasma. The calculation for the partition 
function U is done in three steps. The first is an actual summation over all 
levels belonging to principal quantum numbers for which the energy levels are 
completely tabulated; the second step is a summation through principal 
quantum number 15, assuming hydrogenic levels. In the last step the remaining 
sum is approximated by an integral and all levels x^ are replaced by the 

ionization potential. The contributions of parent configurations other than the 
series limit are included and the final equation is written in a form that 
explicitly shows the effect of higher parent terms in the calculations. 

Numerical results are presented for four different species at several 
temperatures and pressures; these are compared with calculations of other 
authors. In the case of A1 II (Mg I isoeleetronic sequence) , which has two 
parent configurations, the error incurred by considering only the main parent 
term in the calculations is clearly demonstrated. 



INTRODUCTION 


To fully interpret the experimental data obtained from a magnetograph 
system, the theory for the formation of a Zeeman-sensitive absorption line 
in a magnetic field must be investigated. Moe [ 1] has shown how the radiative 
transfer equations for polarized light, as derived by Unno [2] , can be numer- 
ically integrated using any desired atmospheric model. In his derivation Moe 
has made the assumption of local thermodynamic equilibrium so that the rate 
equations and transfer equations are uncoupled and the level populations are 
governed by the Saha -Boltzmann equations. These latter equations require the 
calculation of internal energy partition functions for atoms and ions, taking into 
consideration their dependence on the local temperature and electron pressure 
of the surrounding plasma. This report outlines the program that was developed 
in order that partition functions could be calculated for any atom or ion, using 
any desired atmospheric model. 

PARTITION FUNCTIONS OF ATOMS AND IONS 
IN A PLASMA 

Partition functions enter into astrophysical calculations in the deter- 
mination of atomic level populations through the Saha and Boltzmann equa- 
tions. In particular, one must calculate the internal energy partition function 
for atoms or ions with all energies measured relative to the ground state: 

00 - X-/kT 

U = £ g e 1 . (1) 

i=i 

Here g^ is the statistical weight of the ith level, x^ is the excitation 

energy of the ith level, T is the absolute temperature, and k is the Boltzmann 
constant. For an isolated atom or ion, the partition function is divergent 
because of the combining factors of an infinite number of levels and a series 
limit. However, when an atom (or ion) is immersed in a plasma, the higher 
levels that were bound states in the isolated atom become continuum states 
because of the interaction between the atom and neighboring ions. Thus the 
number of bound levels is finite and the partition function sum is convergent. 
This concept can be analytically expressed as a probability function, P , 

which denotes the probability that an atom in the ith energy level is not 
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perturbed by the surrounding plasma. For the lowest-lying levels, P. is one; 

for the higher levels, the energy of the states decreases, and the optical 
electrons are more influenced by perturbations from surrounding ions on the 
potential field in which they move so that P decreases. Finally, when the 

interaction energy is comparable to the ionization energy of the bound optical 
electrons, the state is destroyed and P^ is zero. The partition function can 

thus be written 


00 - x./kT 

u = z gi v 1 

i=i 


(2) 


The determination of the probability P^ has been approached in several 

ways. The method outlined below is due to Unsold [3] . A different approach 
to the problem is to consider the reduction of the ionization limit as outlined by 
Griem [4] . Basically, Unsold considers P to be given by the probability 

that the nearest neighboring ion is beyond a critical distance R : 

c 

- ~ N R 3 

T-. _ 3 O C / n \ 


where N is the number density of perturbing ions. For singly ionized 
o 

perturbers, N equals the number density of free electrons N ; in general 

O 6 

for a perturbing ion of charge Z ' , 


Z’N = N = P /kT , (4) 

o e e 

where P g is the electron pressure. 

Following the discussion of Pecker and Schatzman [5] , R is that 

separation of the perturbing ion and atom at which the maximum potential 
energy of the optical electron in the resulting field of the atom and ion equals 
the ionization energy of the electron in its orbit. 

From Figure i, with R and r fixed, one can see that the potential 
energy of the electron will be greatest along the line of separation R between 
the atom and perturber. Thus, 
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(5) 



where p+r=R. The effective charge 
of the atom minus the optical electron 
is denoted by Z . If we maximize 

Ep under the constraint that 

p + r = R, with R fixed, there re- 
sults: 


Figure i. Schematic diagram showing 
the relative separations of atom, 
ion and electron. 



( 6 ) 


where r t and p t are the respective 
values of r and p at which Ep is a 


maxi mum . T hus , 


■I 2 


E„ 


P /maximum 


r i 


Pi 


(7) 



Then R^ is the separation distance at which this maximum potential 

energy equals the ionization energy of the optical electron with shell number 

n : 
i 


Z e 2 _ Z'e 2 
r l ~ Pi 


hcRyZ 2 
2 


n 


( 8 ) 


In this expression hydrogenic levels have been assumed for the optical electron; 

Ry is the Rydberg constant, h is Planck's constant and c is the velocity of 

light. Using the relation R =v i +p l and equation (6) , the expression for 

c 

R becomes 
c 


R = n. 2 
c 1 


e 2 ( Z + Z' + 2 n/z Z') 
h c" Ry Z 2 


(9) 
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The probability function P is now written 


P.= e 

1 


C (Z) P 8 n 6 
e 1 


( 10 ) 


where 8 = 


5040 


. The constants C(Z) are tabulated in Alter [6] 


REDUCTION OF THE PARTITION FUNCTION EQUATION FOR 
PRACTICAL CALCULATIONS 


For an atom (or ion) , the partition function is to be calculated from 
the expression 

°° -C (Z) P 6 n. 6 - x./kT 

U = Z g,e 6 1 e 1 . (11) 

1=1 

The energy states y. are tabulated in tables of atomic energy levels [ 7] , but, 

for most atoms, these tables are incomplete, with only levels of the lower 
orbital angular momentum configurations being known. Following the approach 
outlined by Aller, Elste and Jugaku [8] , U is written as the sum of two 
parts: 


U = U 0 + Ui 


(12) 


If the atomic energy level tables become incomplete at principal quantum 

number n. = m , 
x 


U ft 


n.=(m-i) 

= 1 E 

i=i 


g i P i 6 


*i 


/kT 


(13) 


with g. = 2 J + 1 . In the tables of Reference 7, the energy levels are 
expressed as wave numbers so it is convenient to write 

. _4 ix ~ 

y./kT = (2.8548 x 10 ) 9 v\ scj 8 . The term Uj is the remaining 

summation over higher levels of which wave numbers are incompletely 
tabulated. To include all the states of these higher levels, it seems best to 
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assume hydrogenic structure for the levels and perform the sum, rather than 
use the tabulated levels and risk omission of states with high statistical weights. 
Accordingly, 


Ui 


z 

n.=m 


-Cj 9 

'iV 




(14) 


where Ry = 109 678.758 cm 1 , Z = i for neutral atoms, Z = 2 for singly 

rsJ 

ionized atoms, etc. , and p is the wave number limit of the atom or ion. The 

JLi 

statistical weight g. of a hydrogen-like system consisting of a ’’running" 
electron with outer shell number n added to a "parent configuration" with a 
multiplicity of states (2S+ 1) (2L+ 1) is given by 


g = 2 n. 2 (2S+ i) (2L+ i) . 


(15) 


The quantum numbers S and L refer to the parent configuration of the atom 
(or ion) . Usually the term "parent configuration" refers to the ground state of 
the next higher ionization stage; however, some systems have more than one 
parent term and this must be considered in the summation of U^. Suppose there 
are N parent terms listed in the tables of predicted terms . Then 


Ui 


Z 


j=i 


z 

n.=m 


" c i 

(g i»j p i e 


(v L ). - Ry 



(16) 


The wave number limits of the higher parent terms can be written in terms 

of the series limit v : 

Li 


‘Vj ■ 


v + (Ar). 

1 


r*j 


where ( Av) refers to the listed wave number of the 

in the tables for the next higher ionization stage. The 
(g|) j can be written as 


jth parent configuration 
statistical weights 


6 



Then equation ( 16) becomes 


tb = 


x 


2(2 + 1) (2 L t + 1) Yj P i e 


Ci d[ v L - R y ^r 


n =m 

l 


i* Z 


N (2 S. + 1) (2 L. + 1) - Ci 0 (A v ) . 


J 


1 


j=2 


(2 S t + 1) (2 Lj+ 1) 


(17) 


The last term in brackets in equation ( 17) represents the contribution of all 
parent terms other than the series limit term (with quantum numbers Si and 
Li) . This contribution is independent of P and can be easily evaluated for 

a given range of 0 to determine the importance of higher parent terms in the 
calculation of Ui. 


The summation over the hydrogenic levels in Ui is performed through 
n. = 15 (Uii) ; the remaining sum (U 12 ) is approximated by an integral where 

the following assumption is also made: 


v -r - Ry — T ^ , for n. > 15. 

L n. L ’ i 


Thus the expression for U 12 can be written as 


- v 0 °° - [C (Z)P 0 ] n* 

U 12 = 2 (2 Si + 1) (2 Li+ 1) 10 f n 2 e e 

n=i6 


dn 


( 18 ) 


where x T is expressed in eV. With the substitution 


n 


= y 1 / 3 [C(Z) P 0] 1//<3 


there results 



U„ - + g (2 ^ fe-^dy 

3^C (Z) P 0 y 0 

t/ 


(19) 


where y 0 = 16 3 \l C (Z) P 0* . The final result for the total partition 

0 

function becomes: 

2(2 Si + 1) (2 L t + 1) 


U 


= E ( 2 J- + 1) P. e 
1 1 
i=l 


Ci 9 v. 


+ 


x 


15 



n =m 

l 




+ 



3 


4 C (Z) P 9 


(2 Si+ 1) 


" X T 


x (2 Li + 1) 10 


Erfc (y 0 ) 


N (2 S. + 1) (2 L. + 1) 
+ p 2 (2Si+ 1) ( 2 L 1+ 1) 


Ci 0(Av) . 


( 20 ) 


- C (Z) P 0 n. G 

c i 

with P. = e . The error function can be evaluated from the 

i 

expansion 

E rfc (y 0 ) = e y ° (a t t + a 2 t 2 + a 3 t 3 + a 4 t 4 + a 5 t 5 ) , 

where t = [1 + (0.3275911) y 0 ] 1 , and 

a t = 0.254829592, 
a 2 - -0.284496736, 
a 3 = 1.421413741, 
a 4 = -1.453152027, 
a 5 = 1.061405429. 
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NUMERICAL CALCULATIONS 


A FORTRAN program has been written specifically for the calculation 
of Fe I and Fe II partition functions, which occur in the analysis of the forma- 
tion of the Fe I 5250. 2-A line in a magnetic field, following the approach of 
Moe [ 1] . The program was designed to be flexible enough to be adapted to 
other atoms or ions, as well as to be able to handle various photospheric and 
sunspot models . The program is divided into three parts corresponding to 
the calculations of U 0 , Uji and U 12 . 

In the calculation of U 0 , all the wave numbers ( r.) , statistical 

weights (2 J.+ i) , and shell numbers (n.) must be read into the program. 

The shell number m at which the wave number tables become incomplete 
can be determined by a comparison of the tables of the observed and predicted 
terms. If the observed terms are almost complete for a given shell number, 
the missing states can be estimated from corresponding terms of other mem- 
bers of the isoelectronic sequence. Or, if the missing states have small 
statistical weights and correspond to higher wave numbers, they may be 
neglected in the sum. The calculations for 1% and U 12 are straightforward. 

To check out the program, partition functions were calculated for four 
different species, ranging from simple-series spectra to complex spectra: 

A1 IH (Na I isoelectronic sequence) , A1 II (Mg I sequence, two parent con- 
figurations) , Fe I (Fe I sequence, 19 parent configurations) , and Fe II 
(Mnl sequence, 24 parent configurations) . 

Al III 


For Al HI, the tables of wave numbers are complete from the 3s 2 S 
through the 6h 2 H° configuration; the sum in is thus started at n = 7. 

In compiling the known levels, each substate was considered separately for 
the three lowest levels. The parent configuration of Al III is a *S 0 state, 
which indicates a parent multiplicity of one. The partition functions for Al III 
were calculated for various electron pressures at die temperature 9 - 0. 14; 
the results are listed in Table 1 together with Al III partition functions cal- 
culated by Aller, Elste and Jugaku [6,8]. 
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TABLE 1. COMPUTED PARTITION FUNCTIONS FOR AI III 


e 

Log P^ 

AI HI U 

AI III U 3 

0. 14 

2.0 

13.086 

13. 18 

0. 14 

3.0 

6. 16 

6. 17 

0. 14 

4.0 

3.97 

3.98 


a. [6,8]. 


AI II 

For singly ionized aluminum, AI II, the calculations are more com- 
plex because of the presence of two parent configurations in the table of 

observed terms: 2 S and 2 P° . For the n. = 3 shell number of the optical 

i 

electron, there are four missing terms: the *S term of the 3p 2 configuration 
and the 1 P° , 1 D° and 1 F° terms of the 3p ( 2 P°) 3d' configuration. Wave 
numbers were estimated for these levels using the wave numbers of corre- 
sponding terms of Si III and P IV. 

Because of the existence of the two parent terms, the partition function 
was evaluated twice, first considering only contributions for n. ^ 4 from the 

levels of terms belonging to the 2 S parent configuration, and then considering 
contributions of both parent configurations. For the first case of the 2 S parent 
term alone, the observed terms are complete through shell number five so 
that the sum U n was begun at n. = 6. The substates of the first seven terms 

of U 0 were summed separately; thereafter the terms were grouped as single 
levels. The multiplicity (2S+ i) (2L+ 1) for the series limit 2 S is 2. 

The resulting values for the partition function calculated in this manner are 
listed in Table 2 as U ( 2 S) . 


TABLE 2. COMPUTED PARTITION FUNCTIONS FOR AI II 


0 

Log P 

e 

U( 2 S) 

U( 2 S, 2 P°) 

u a 

0. 14 

2.0 

213.2 

285. 7 

213.8 

0. 14 

3.0 

72.2 

95.6 

72.4 

0. 14 

4.0 

26.3 

33.7 

26.9 


a. [6,8]. 
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When the contributions of both parent terms were considered, the 
sum in U u was initiated at n. = 4 since the observed terms of the 2 P° 

i 

parent configurations are very incomplete for n. > 3. The partition function 
thus takes the form 


- c 1 6 v. 

U = Yj ( 2 J . + i) p - e 


n=3 

l 


-c^b-Ry-2 

x P e \ L i /+ 
i 


15 


2(2 St + 1) (2 Lj+ 1) 2 


n. 


n =4 

l 


n/T 

3 


si c (z) p e 

e 


Erfc (y 0 ) (2SJ+1) (2L t +l) 10 


V 1 


x 


For A1 II, 


( 2 S? + 1) (2 L9 + 1) - Ci 6 A i >2 
1 + (2S t + 1) (2 Li+ 1) 6 

— 53 860 cm 1 so that 


( 21 ) 


(2S, + 1) (2 L 9 + 1) - c x e Ar 2 - 15.376 9 

(2 Si + 1) (2 Li + 1) 6 6 ; for 9 = 0. 14, this 

contributes a factor of 0.35 more to [1% + U 12 ] . Numerical results using 
equation (21) are presented in Table 2 under the heading U( 2 S, 2 P°) . The 
last column represents data published by Aller, Elste and Jugaku [6,8], The 
effect of contributions from other parent terms can be clearly seen. 

Fe I 


The Fe I and Fe II partition function programs were initiated to 
calculate the populations of the lower level of the Fe I 5250.2 A transition 
for various photospheric and sunspot models for which 0.4 < 9 < 1. 9. 

Thus the effects of extremely high temperatures have been ignored. For Fe I, 
the sum U 0 included the 3d 6 4s 2 and 3d 8 levels plus all listed terms for 
which the active electron had a shell number n. = 4. Missing states in the 

o 

3d configuration were estimated from similar terms in Co n and Ni III 
and from similar terms in members of the Ca I isoelectronic sequence, 
assuming a similarity between die 3d 8 and 3d 2 configurations. The missing 
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3d 6 4s 2 states were estimated from known wave numbers of tabulated members 
of that configuration. Missing levels for the m = 4 shell number were 

numerous and no attempt was made to estimate them; rather, for a represent- 
ative 6 and P , their total contribution to the partition function was estimated, 

6 -4 

using hydrogenic levels, and found to be « i.O x 10 . Therefore, no cor- 

rections were deemed necessary for them and the sum was initiated at 
n. = 5. For both and U 12 , only the series limit configuration a e D was 

considered as a parent term, although 19 actual parent terms exist. The 
effect of omitting these other parent configurations can be determined from the 
last term in brackets in equation ( 17) . For a representative temperature for 
the solar photosphere, 0= 0. 9705, this term has the value ( 1 + 0.63) . How- 
ever, for the a 6 D parent term alone (Uu + U^) = 0.016 (assuming P^ = 

3.554 dynes/cm 2 ) . Thus all parent terms contribute a value of 0.026 to the 
total partition function. But at this temperature and pressure, U 0 = 28.486, 
so that it seems justified to consider only the a e D parent term for photospherie 
and sunspot models; the contribution of the other 18 parent terms is 
insignificant, percentage wise. Table 3 compares the results of the computer 
program with tabulated values of Muller and Mutschlecner [9] . 

Fe II 

From the tables of observed and predicted terms for Fe II, it can be 
seen that the 3d 7 configuration is complete. Most of the remaining tabulated 
levels belong to 3d 5 4s 2 or 3d 6 ( parent) 4x configurations (x = s, p, d, f ) , 
which are incomplete. Numerical estimates for 9 = 1.00 and 
P g = 2.5 dynes/cm 2 indicate that all missing levels of the 3d 5 4s 2 configura- 
tion would contribute « 2. 5 x 10 3 to the partition function. Similar estimates 
for all other missing states with outer shell number n. = 4 indicate their 

—3 

total contribution would be less than 1.3 x 10 . Consequently the 1% 

sum was begun at n. = 5 and the U 0 sum included all listed levels with outer 

shell number n. £ 4. In the calculations for U 1; t and U 12 , only the main 

parent term a 5 D was used; contributions of the other parent configurations 
are small at the temperatures considered because of the large wave number 
gap between the a 5 D ground-state term (v 1027 cm ) and the next lowest 
term, a S P (~ « 20 000 cm ), in the Fe III level structure. The term in 
equation ( 17) that represents the contributions of other parent terms to 
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TABLE 3 . PARTITION FUNCTION OF NEUTRAL IRON AS A FUNCTION OF 
TEMPERATURE AND ELECTRON PRESSURE 


Log T 0 

e a 

P 

e 

U 

b 

u 

-3.6 

1.0944 

0. 1323 

26.465 

26.430 

-3.4 

1.0938 

0. 1647 

26.472 

26.438 

-3.2 

1.0933 

0.2044 

26.478 

26.444 

-2.8 

1.0913 

0.3150 

26.505 

26.471 

-2.4 

1.0864 

0.4882 

26.572 

26.539 

-2.0 

1.0748 

0.7733 

26.736 

26.703 

-1.6 

1.0480 

1.3070 

27. 138 

27. 103 

-1.2 

0.9990 

2.4800 

27.961 

27.921 

-1.0 

0.9705 

3.5540 

28.502 

28.457 

-0.8 

0.9400 

5.3040 

29. 140 

29.089 

-0.4 

0.8703 

14.790 

30.888 

30.813 

0 

0.7835 

62.450 

33.868 

33.723 

0.4 

0.6845 

356.30 

39. 129 

38.760 

0.8 

0.5688 

2788.0 

50.538 

49. 195 

1.2 

0.4365 

26 780 

82. 113 

75.980 


a. 6 = 5040/ T 

b. [9], 
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sum 11^+ U 12 can easily be shown to be ~ 0. 01 for 0= 0. 9705 and 

P^ = 3 . 554 dynes/ cm 2 ; the a 5 D term itself contributes practically nothing 

to the partition function at this 0 and P so, again, the neglect of these 

6 

other terms is justified. Sample calculations from the computer program are 
given in Table 4 along with values from the paper of Miiller and Mutschlecner 
[9]. 


CONCLUSION 

In conclusion, it should be emphasized that for higher temperatures 
the contributions of all parent terms should be carefully examined. In a 
paper by Traving, Baschek and Holweger [ 10] , these effects have been con- 
sidered in the calculations of partition functions for atoms and ions of 
astrophysical interest. The resulting partition functions are presented in 
polynomial form for easy machine calculations, eliminating the necessity of 
reading the tabulated energy levels into the program. 
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